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We review and compare different variational formulations for the Schrodinger field. Some of 
them rely on the addition of a conveniently chosen total time derivative to the hermitic Lagrangian. 
Alternatively, the Dirac-Bergmann algorithm yields the Schrodinger equation first as a consistency 
condition in the full phase space, second as canonical equation in the reduced phase space. The 
two methods lead to the same (reduced) Hamiltonian. As a third possibility, the Faddeev-Jackiw 
method is shown to be a shortcut of the Dirac method. By implementing the quantization scheme 
for systems with second class constraints, inconsistencies of previous treatments are eliminated. 
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Outstanding equation of modern physics, the Schrodinger equation has multiple and deep connections with integral 
principles. Historically, Schrodinger obtained his equation guided by the beautiful analogy between the Fermat 
principle and the principle of least action [1]. Moreover, motivated by a remark of Dirac [2], Feynman has derived 
the Schrodinger equation from the Huygens principle, realizing the first step towards his path integral approach [3], 
[4]. (For recent recent reviews see [5], [6].) The Schrodinger field is equally a popular choice to illustrate how second 
quantization proceeds. 

| At a closer inspection, however, the power and beauty of the variational approach is obstructed by an aesthetical 
■ bug. The reason: the abundance of dynamical variables, some of them being redundant. The Lagrangian yielding 
the nonrelativistic Schrodinger equation is linear in the time derivatives of the fields \& and its complex conjugate , F* . 
Thus the Legendre transformation does not lead to an unambiguous Hamiltonian. Various Hamiltonians, all having 
Schrodinger's equation as canonical equation can be found in textbooks. They depend either on two pairs of canonical 
variables [7], or - as a result of adding a total time derivative to the hermitic Lagrangian - just on a single (complex or 
real) canonical pair [8]- [10]. In the next section we briefly rewiev these approaches, pointing out both the weakness 
| and ingenuity in sweeping away the problem. 

There are two equivalent methods to circumvent this difficulty. In section 3 we treat the Schrodinger field as a 
constrained system, applying the Dirac-Bergmann algorithm [11]- [13]. As the system has second class constraints, 
JL . the dynamics involves Dirac brackets. We present two alternatives to the existing derivations of the Schrodinger 
equation. First, the consistency requirement yields the Schrodinger equation in the form of a weak equation on the 
full phase space. Second, by a suitable canonical transformation we introduce new canonical coordinates, containing 
the constraints. The Dirac bracket of the full phase space becomes the Poisson bracket of the reduced phase space 
and the Schrodinger equation is a canonical equation. As a bonus we recover the Hamiltonian obtained by addition 
of a properly chosen time-derivative. At the end of Section 3 we present an alternative discussion in terms of real 
variables: the real and imaginary parts of the complex field 

In Section 4 we apply the Faddeev-Jackiw scheme developed for Lagrangians, which are first order in "velocities" 
[14], [15]. We verify that the fundamental brackets of the Faddeev-Jackiw approach coincide with the Dirac brackets. 

Finally in the fifth section we follow the canonical quantization scheme [11]- [13], [16], giving an operator repre- 
sentation of the Dirac bracket algebra of the canonical variables. This is equivalent with the quantization of the 
Faddeev-Jackiw fundamental bracket. Our approach avoids the interpretational inconsistencies of some standard 
treatments [7], [8], already pointed out by Tassie [17] and is a viable alternative to the existing reduced phase space 
quantization schemes [9], [10]. By imposing the second class constraints as operator identities, second quantization 
proceeds smoothly. 

II. STANDARD VARIATIONAL PROCEDURES. A REVIEW 
A. Two pairs of complex variables 

The action for the Schrodinger field cf. Henley and Thirring [7] is: 

S[ip,?p*} = J dtj dr C (2.1) 
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£ = V (ip*ip-ipTp*) - —Vip*Vip-V(r,t)ip*ip • 
2 V / 2m 

This Lagrangian density C is hermitian 1 . Variation of (2.1) with respect to ip* and ip gives the Schrodinger equation 

• h 2 

ity + ^ A V> -Vil> = , (2.2) 
2m 

and its complex conjugate. The canonical momenta of the complex conjugate variables are complex conjugates too: 

7T = f r , 7T* = • (2-3) 

According to Henley and Thirring, the Hamiltonian density H, when written in terms of the phase space variables 
(ip, ip*, ^ t*)) to be: 



S[ip,Tr,ip*,TT*\ = f dt f dr (W> + T*f - W) 



W = ^-(VV*Vtt* - V^Vtt) - Jy(r,t)(^7r- V*tt*) • (2.4) 

Variations with respect to "0, 0>* give the Schrodinger equation and its complex conjugate. The same equations 
emerge as a result of the relations (2.3) and the variations with respect to n, ir* . The Schrodinger equation as such 
a frequent outcome indicates that there are redundant variables in the formalism. 

Note that it is impossible to express the velocities ip, ip* m terms of canonical data from the expressions of momenta 

(2.3) , as would be required by the Legendre transformation. How can one then find the Hamiltonian (2.4)? No hint for 
this is given in [7]. We can still invent a method. First express the fields ip, ip* from (2.3), insert them in the kinetic 
terms of (2.1), but not in the time derivatives. Now the action (2.1) takes an already Hamiltonian form, similarly as 

(2.4) . At this point we see a Hamiltonian expressed in terms of canonical data, which is hermitian, however does not 
give the correct equations. Instead we devide the potential terms of (2.1) in two equal parts, eliminate the starred 
fields from one part and the unstarred ones from the other part by means of the relations (2.3). The hermiticity of 
the Hamiltonian is kept and we get the Schrodinger equation. Proceeding in other ways, for example, eliminating 
all starred or all unstarred fields (and consequently destroying the hermiticity) will result in wrong Hamiltonians. 
However, as we have seen, hermiticity alone is not a criteria for correctness. 



B. One pair of complex variables 

Furthermore, the requirement of hermiticity is not even a necessary one. Schiff [8] derives the Schrodinger equation 
from a non-hermitic Lagrangian, found from (2.1) by adding the total time derivative in(tpip*) '/2: 



S[ip,ip*] = J dt J dr C s 
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C s = ihip*ip-—Vip*V*IJ-V(r,tU*il>. (2.5) 
2m 

Variation of the action (2.5) with respect to ip and ip* gives the Schrodinger equation (2.2) and its complex conjugate. 
However the canonical momenta are not complex conjugate any more: 

tts := —j- = ihip , ttj := -j- = . (2.6) 

dip dip 

The first of these equations is used to eliminate ip* from the action, which becomes: 



1 The Lagrangian density (2.1) when written in terms of the modulus and argument of ip, leads to equations very similar to 
the equations of hydrodynamical flow in presence of a potential V . This analogy stands at the base of the hydrodynamical 
model of quantum mechanics. 
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S[iP,n s ] = fdtjdr (tt s ^-H s ) 



H s = -^VV'Vtts- l -V{r,t)i>K S • (2.7) 
2m n 

Thus no starred field shows up as canonical variable. The canonical equation obtained by varying ip is the Schrodingcr 
equation. The other canonical equation from variation of its, together with the first relation (2.6) gives the complex 
conjugate Scrodinger equation. Some of the superfluous variables were removed by the addition of a total time 
derivative, achieving a partial reduction to the true degrees of freedom. This Hamiltonian description is the simplest 
one in terms of complex fields. 



C. One pair of real variables 



A more efficient way to look on the variational problem for quantum mechanics is described by Kuchaf [9] and 
alternatively by Cohen- Tannoudji, Dupont-Roc and Grynberg [10]. These approaches rely on the decomposition of 
the field tp in real and imaginary parts: 



?= -^(V' + V'*) , P = 
In terms of the real fields q and p the action (2.1) becomes: 



(2.8 



S[q,p] 



dt dr C 



(Vg) : 



Vp 
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(2.9) 



By adding the total time derivative (pq)' /2 to the Lagrangian density, the action takes an already Hamiltonian form. 
The field p turns to be the conjugate momentum to q: 



S[q,p] 



dt 



Am 



dr (pq - Hk) 
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The two canonical equations obtained by variations with respect to q and p are: 



1 a V 
2m h 



P= 7T- Aq-Vq . 
2m 



(2.10) 



(2.11) 



Up to global factors, the first equation of (2.11) is the real part, while the second the imaginary part of the Schrodingcr 
equation (2.2). 

In this latest approach a full reduction of the phase space to the true degrees of freedom was achieved, as the 
Schrodinger equation emerges only once in the Hamiltonian formalism. In terms of real fields this is the simplest 
description, which again relies on the addition of a properly chosen total time derivative to the Lagrangian density. 

We will see in the next section that the addition of specific total time derivative terms to the Lagrangian (2.1) is 
not compulsory. The standard Dirac-Bergmann algorithm leads directly either to the Hamiltonian density (2.7), (in 
a description in terms of complex fields) or the Hamiltonian density (2.10) (if real fields are introduced). 



III. THE DIRAC-BERGMANN ALGORITHM 



Constrained systems are characterized by the singularity of the inertia matrix, whoose elements are given by the 
second derivatives of the Lagrangian with respect to the generalized velocities. This property always holds if a 
Lagrangian density is linear in the time derivatives of fields [13], as is the case for (2.1). Then the Dirac-Bergmann 
algorithm takes the role of the Legendre transformation. 
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A. Complex fields 



The momenta (2.3) provide two primary Hamiltonian constraints: 

<Ai := Tr-yV* = 0, 2 := 7r * + ^ = O. (3.1) 

Time evolution of an arbitrary phase space function / is generated through the Poisson bracket by the primary 
Hamiltonian density Hp rather than the canonical Hamiltonian density He ■ 



dr Hp 



f *{f,H P } , H P = J 
Uv ■= Hc+ip<pi+ip*<p 2 

h 2 

He := — VV'*VV' + ^(r,i)V*V • (3.2) 

777 

Here the velocities ip, ip* are unknown functions of the phase space variables. The symbol S3 denotes weak equality, 
holding only on the hypersurface determined by the constraints (3.1). 

Consistency requires the time derivatives of the primary constraints <pi, 4>2 to vanish: 

h 2 

S3 (Pi w + AV>* - vv* 

2m 

?? 2 

S3 02 w ifitp + - — Atp - Vip* . (3.3) 
2m 

These are the Schrodinger equation and its complex conjugate. They emerge as weak equalities; a price one has to 
pay for working on the complete phase space (tp, ir, ip*, n*). No secondary constraint appears in the theory since 
Eqs. (3.3) are relations determining the unknown functions ip and tp* . 

The Poisson bracket of the two constraints shows that they are of second class: 

{01,0a} :={0iM),0 2 (r',i)} = -iM(r-r'). (3.4) 
Time evolution can be given in the alternative form: 

./ = {f,H c } D , H c = J dvHc (3.5) 

in terms of the Dirac bracket [11]: 

U,9}d := {f,g}- Yl {f,<t>i)({<l>k,<t>l})^{<l>j,9} ■ (3-6) 

i ,3=1,2 

Here ({</>&, <Pi\) 1 denotes the inverse of the matrix with elements given by the Poisson brackets of the constraints. 
Straightforward computation using (3.4) gives the following expression for the Dirac bracket: 

u.,u - I{/„ } -i/*(g* - g + ( 3.r) 



From here it is immediate to write the Dirac brackets of the canonical data: 

i 

n 



{tP(r),r(r')} D = -U(r-r') , {tP(r),iP(r>)} D = {V>*(r), PV)}d = , (3i 



if) 

{n(r),n*(r')} D = -jJ(r - r') , {n(r), 7r(r')} B - K(r), 7r*(r')}i> - 0, (3.9) 
{tP(r),n*(v')} D = {r(r),n(r')} D = 0, (3.10) 
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Wr),n(r')} D = \{^{v), 7r(r')} = \s(t - r') , 



(3.11) 



{^(r),7r*(r')} D = i{V*(r),7r*(r')} - ^(r - r') . (3.12) 

These Dirac brackets do not contain phase-space functions, thus no operator ordering difficulties will occur during 
quantization. 

Reduction. Dirac brackets of second class constraints with arbitrary functions vanish. Thus the constraints can 
be solved prior to calculating the Dirac brackets, by reducing the phase space to the physical degrees of freedom. As 
each second class constraint reduces the dimension of the phase space by one, a basis in the reduced phase space is 
provided by a single pair of complex canonical data. A suitable canonical transformation turns the constraints into 
the other pair of canonical data: 

/ V V \ _^ ( *l>i = ip/2+in*/h V2 = -i<h/K 

y-K TT* J y 7Ti = 7T + ihtp* /2 7T 2 = 01 

A straightforward check shows that the Dirac bracket (3.7) written in terms of the new coordinates becomes the 
Poisson bracket of the reduced phase space, coordinatized by tpi , 7Ti : 

{Lg}D = J dr {Wi5^~s^Wj • (3 ' 14) 

Now it is immediate to verify the generic property that the Dirac brackets of the constraints V>2 and 7r 2 with arbitrary 
functions vanish. 

The Hamiltonian density on the reduced phase space is the one introduced by Schiff (2.7), with if>i, tt\ in place of 
ip and tts- The canonical equations are the Schrodingcr equation and its complex conjugate. 



(3.13) 



B. Real fields 



Our starting point in this section is the action (2.9) written in terms of the real fields q and p. As in the previous 
section we have presented in detail the method, here we merely list the results. From the definition of the momenta 
P q . p canonically conjugate to the q,p variables we find the constraints which are second class: 



$1 == P q ~ f , $2 := P P + § , 
{$ 2 ,$!}:= S(r-r'). (3.15) 



The consistency requirements $1.2 ~ are the real and imaginary parts of the Schrodinger equation, Eqs. (2.11). 

Reduction. The constraints (3.15) already form a canonical pair of variables, thus a canonical transformation to 
the reduced phase space is immediate: 



q P \ ( Qi = 9/2 - P P Q2 = $2 

P g P P \Pi=p/2 + P q P 2 = $i 



(3.16) 



Again the Dirac bracket on the full phase space (Qi^, Pi, 2) becomes the Poisson bracket on the reduced phase space 
(<5i, Pi). The Hamiltonian density on this reduced phase space is Hie given in (2.10), with Qi,Pi in place of q,p. So 
the canonical equations on the reduced phase space are the real and imaginary parts of the Schrodingcr equation. 



IV. THE FADDEEV-JACKIW APPROACH 



Developed as a Hamiltonian formulation of dynamical systems with Lagrangians linear in velocities, the Faddeev- 
Jackiw method provides undoubtedly the shortest path toward a fundamental bracket in the phase space. 

The kinetic part of the Lagrangian (2.1) with 2x 00 basic variables C( r ) = (V^ 1 "')) V'*( r ")) determines the symplectic 
2-form with the inverse 

^ (r ' r ] ~ ih \-S(r'-r") ) ' 
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in terms of which the time evolution of the fundamental variables is 

?W=» ii MstfP) H c- (4-2) 

This is nothing but the Schrodinger equation and its complex conjugate. Eq. (4.2) represents a Hamiltonian evolution 
if the fundamental bracket obeys 

{f(r),£V)}FJ = w tf (r,r') , (4.3) 
which is nothing but a shorthand notation for Eqs. (3.8). 

V. SECOND QUANTIZATION 

The standard procedure for canonical quantization of systems with second class constraints is to turn the Dirac 
brackets into commutators cf. the scheme: 

{L 9)0 = 1 -> [f,g} = ihi. (5.1) 

Here /, g and I are phase-space functions, /, g and I are operators. The complex conjugate /* of a function / becomes 
the adjoint operator p. When we apply these prescriptions to the Dirac brackets (3.8)-(3.12), we get: 

fo£(r),^(r')] = S(r-r'), $(r),^(r')] - [^(r),^( r ')] = 0, (5.2) 
[7r(r),7rt(r')] = - ^(r - r') , [7r(r), 7r(r')] = {^(r), ttV)] = , (5.3) 

[V<r),7r(r')] - [^(r), ^(r')] - ^(r-r'), $(r),7rV)] = [^W^(r')] - 0. (5.4) 
The second class constraints (3.1) of the theory become operator identities: 

7T=yV, = • (5-5) 

By inserting Eqs. (5.5) in Eqs. (5.3) and (5.4) we recover again Eqs. (5.2). 

The commutators (5.2) between ip and tfj^ represent the starting point in the second quantization, as indicated by 
Henley and Thirring [7] and described in detail by Schiff [8]. While Eqs. (5.2) have emerged in a natural way from the 
Dirac bracket quantization, they had to be imposed "by hand" in the previous approaches, a feature already criticized 
by Tassie [17]. Schiff arrives to Eqs. (5.2) by imposing the canonical commutation rules: 

[VU S ] = [$+,4] = iM(r-r'), (5.6) 

however his treatment also requires = as can be seen from Eq. (2.6), which is in obvious contradiction with 
the second commutator (5.6). Meanwhile, the treatment of Henley and Thirring starts by postulating Eqs. (5.4), 
in other words with the surprizing statement that what was canonically conjugate in the classical theory in not any 
more canonically conjugate in the quantum theory (see the extra factor of 1/2). They impose the commutators (5.4) 
motivated by the analogy with the variational problem of the harmonic oscillator, written in complex coordinates. 
Tassie proposes a solution to these conceptual problems. By working on the momentum space, he esentially eliminates 
the imaginary part of ip and gives a description in terms of real fields without encountering the above-mentioned 
inconvenience. 

No problems appear in the approaches employing real fields. This is because any description in terms of real fields 
[9,10,17] esentially means that we have eliminated the redundant variables, thus we are quantizing on the reduced 
phase-space. 

No difficulties appear in the Dirac bracket quantization either. The Dirac bracket (3.11) of the variables ip and 
7r, canonically conjugate at the classical level, is one-half their Poisson bracket, thus no reason to wonder why the 
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corresponding commutator contains the factor of 1/2. Consistent canonical quantization of the complex Schrodinger 
field requires the Dirac bracket. 

Stated in other way, if we start from the canonical chart (3.13), compute the Dirac bracket cf. Eq. (3.14) and apply 
the prescription (5.1) for the variables spanning the reduced phase space, we find: 

{Vi(r),7ri(r / )} D = S(r-r') -> [^(r), ^(r')] = iM(r-r'). (5.7) 

But modulo the constraints this is consistent with Eq. (5.4). Now, in contrast with the treatment of Schiff, we can 
impose ip2 = ^2 = 0, because in the framework of the constrained systems these canonical variables are second 
class constraints and they should not be turned into canonically conjugate operators [11]. 

VI. CONCLUDING REMARKS 

We have reviewed how the Schrodinger equation can be found from various Hamiltonians, representing different 
stages of reduction. Hopefully our treatment shed light on the many Hamiltonian formulations of the Schrodinger 
system and their multiple interconnections. We have seen how the construction of these Hamiltonians requires some 
artwork, like the addition of appropriately chosen time derivative terms to the Lagrangian. 

Alternatively the Schrodinger field appears as a computationally simple example for constrained systems. By 
employing the characteristic toolchest, we have found the Hamiltonian and the Schrodinger equation via the Dirac- 
Bergmann algorithm and the consistency requirement, respectively. We have achieved the reduction to the physical 
degrees of freedom by suitable canonical transformations. The canonical equation in the reduced phase space is again 
the Schrodinger equation. 

We have shown how second quantization of the complex Schrodinger field by turning the Dirac bracket (or the 
equivalent fundamental bracket in the Faddeev-Jackiw approach) into commutators on the one side avoids interpre- 
tational difficulties, on the other side leads to the same quantum theory, which emerges from quantization on the 
reduced phase space. 

The Schodinger equation completes the list of famous equations of modern physics, like Maxwell and Einstein 
equations, not covered by the "usual" variational treatments. The exceptions turn out to be rather generic. However, 
there is a major difference: the Schrodinger constraints are second class as opposed to the first class constraints of 
electrodynamics and general relativity. Following the equivalent Faddeev-Jackiw approach, there are no constraints 
at all. This latter approach with uncontested simplicity yields the correct fundamental bracket for the Schrodinger 
system, but the role of the previously found Hamiltonians is revealed only by the Dirac method. 
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